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Central limit theorem for functionals of a generalized 
self-similar Gaussian process 

Daniel Harnett* and David Nnalart^ 


Abstract 

We consider a class of self-similar, continuous Gaussian processes that do not necessarily have 
stationary increments. We prove a version of the Breuer-Major theorem for this class, that is, subject 
to conditions on the covariance function, a generic functional of the process increments converges 
in law to a Gaussian random variable. The proof is based on the Fourth Moment Theorem. We 
give examples of five non-stationary processes that satisfy these conditions. 


1 Introduction 

We consider a centered Gaussian process X = {Xt,t > 0} that is self-similar of order 
(3 e (0,1). That is, the process {a~^Xat-it > 0} has the same distribution as the process 
X for any a > 0. Consider the function 0 : [1, oo) —)■ M given by 

0(a:)=E[XiXj. (1) 

This function characterizes the covariance function. Indeed, for 0 < s < t, we have 

R{s, t)=E [X,W] = 5^% [XiW/d = (;) • 

The best known self-similar Gaussian process is the fractional Brownian motion 
(fBm), where 

and the self-similarity exponent (3 is the Hurst parameter H G (0,1). For the fBm, 

(j){x) ^ ^ — {x — , x > 1. 

We will impose conditions on the function 0 such that E — W)^] ~ s" as s —)■ 0, 

for some constant a G (0, 2(3] that we call the increment exponent. For example, a = 2H 
for the fBm, since E [(W+s — ^t)^] = ■ However, there are examples where a < 2(3. 
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Our goal in this paper is to identify a set of conditions on a, jS and snch that 
we can establish a central limit theorem for functionals of the increments of X. More 
precisely, let 7 = A/'(0, 1 ) and consider a function / G 7 ), which has an expansion 

of the form 

00 

fix) = '^CqHq{x), ( 2 ) 

q=d 

where d > 1, q fy 0, and Hq{x) is the gth Hermite polynomial, dehned as 


Hqix) = 

The index d is called the Hermite rank of /. 
For integers n >2 and j > 0 dehne 


2 d'' 




e 2 . 


AX, 


XX£ — Xi +1 — Xj_ and . ry ■ 

n n n 11 A A j_ 11 ^2 (Q) 


We consider the stochastic process dehned by 


Fn{t) = ^ Y] [nt\ > 1 , 

j=0 


( 3 ) 


and Fn{t) = 0 if \nt\ <1. It is well known that if the process X has stationary 
increments, the convergence to a normal law for the seqnence of random variables Fn{t) 
can be deduced from the following central limit theorem proved by Breuer and Major in 

[3]. 

Theorem 1.1. Suppose {Yj,j > 1} is a centered stationary Gaussian sequence with unit 
variance, and denote by p{k) = E(FkFk+fc) its covariance function. Consider a function 
f G L^(M, 7 ) with Hermite rank d > 1. Then the functional 


n 

= (4) 

j=i 

converges in distribution to a normal law A/'(0, cr^) as n tends to infinity, provided 
EfcezlP(^)l'^ < in this case p(/c)T 

For example, in the case of the fBm, the sequence 0<j<n — 1} dehned by 

- Bf) , 

\ n 71 / 

is stationary, and the Breuer-Major theorem implies that if d > 2 and FT < 1 — the 
sequence of random variables 


n—1 


j=0 





1 


( 5 ) 
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converges in law as n tends to infinity to a Ganssian random variable, with mean zero 
and variance given by 

=i; E (I™+ii“ - 2i“i“+1“ - ii“)’ ■ (6) 

q=d mgZ 

See [3] and Theorem 7.4.1 of [T3] . 

Onr main result (see Theorem 3.4) says that if the covariance of X satishes certain 
conditions and the increment exponent a satishes 0 < a < 2 — 4, then the hnite dimen¬ 
sional distributions of processes {Fn{t),t > 0} dehned in ([3]) converge in law to those 
of a Brownian motion with scaling given by ([6]), where 2H is replaced by a. Notice 
that the sequence of scaled increments {Yj^n,j > 0} is not necessarily stationary and 
we cannot deduce this result form the Brener-Major theorem. On the other hand, the 
relevant parameter in the limit theorem is the increment exponent a, instead of the 
self-similarity parameter /3. 

The convergence in law of the hnite-dimensional distributions in Theorem 3.4 follows 
from the Fourth Moment Theorem [T31 [TB] . which represents a drastic simplification of 
the method of moments to show the convergence to a normal distribution. In order to 
establish convergence, it is sufficient to show the convergence of the variances, and that a 
condition involving the contraction operator is satisfied. In this paper, we use extended 
versions of the Fourth Moment Theorem mm- 

For the particular case of a single Hermite polynomial, that is f = Hq for some q >2, 
one can show the convergence in total variation of the marginal distributions and a 
functional central limit theorem (see Section 3.4). 

We show examples of known processes that satisfy the required conditions, including: 

(a) bifractional Brownian motion (see i), 

(b) subfractional Brownian motion (see j2]), 

(c) an ‘arcsine’ Gaussian process introduced in a paper by Jason Swanson |20], and 

(d) two self-similar Gaussian processes that form the decomposition of a process dis¬ 
cussed in a paper by Durieu and Wang [n|. 

In examples (a), (b) and (d), the self-similarity and increment exponents are the same, 
that is, a = 2/3. This is not true in example (c), where a < 2/3. 

The outline of this paper is as follows. In Section 2, we give definitions and background 
needed to use the Fourth Moment Theorem. In Section 3, we introduce a set of covariance 
conditions on X, and we state Theorem 13.41 which shows that the finite dimensional 
distributions of Fh converge in law when these conditions are met. Some particular 
applications of Theorem 13.41 are discussed. Section 4 discusses the examples (a) - (d) 
above, and Section 5 contains the proofs of two technical lemmas. 

1.1 Acknowledgement 

The authors wish to thank an anonymous referee for helpful comments, especially re¬ 
garding the proof of Theorem 13.41 
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2 Theoretical background 

Following is a brief description of some identities that will be nsed. The reader may refer 
to [131 E] for detailed coverage of this topic. Let Z = {Z{h), h G T-L} be an isonormal 
Gaussian process on a probability space {Q,iF,P), indexed by a real separable Hilbert 
space l-i. That is, Z is a family of Gaussian random variables such that E[Z(/i)] = 0 and 
E [Z{h)Z{gy\ = {h,g)y^ for all h,gE'H. We will always assume that P is the hltration 
generated by Z. 

For integers g > 1, let 1-i®^ denote the gth tensor product of l-i, and let denote 
the subspace of symmetric elements of 

Let {e„,n > 1} be a complete orthormal system in H. For functions f,gE 'H®'^ and 
p G {1,..., q}, we dehne the pth-order contraction of / and g as that element of 
given by 

00 

f®p9= ^ (7) 

i'l 

where f <Zo 9 = f Z) 9 hy dehnition and, if /, g G 'H®'^, f <Zq 9 = (/, In particular, 

if /, g are symmetric functions in /i^) for a measure /i, then we have 

f®i9= [ f{s,ti)9{s,t2) Kds). ( 8 ) 

Jm 

Let Hq be the gth Wiener chaos of Z, that is, the closed linear subspace of L‘^{Q) 
generated by the random variables {Hq{Z{h)),h G "H, ||h||-H = 1}, where Hq{x) is the gth 
Hermite polynomial. It can be shown (see [13], Proposition 2.2.1) that if Z, Y A7(0,l) 
are jointly Gaussian, then 


:x;'. 

(9) 

For g > 1, it is known that the map 



(10) 

provides a linear isometry between 'H®'? (equipped with the modified norm ■\/q\\\ ■ 
and Hq. The random variable Iq{-) is the generalized Wiener-Ito stochastic integral (see 
[T3]. Theorem 2.2.7). By convention, "Ho = M and Io{x) = x. 

It is well known that L^(G) can be decomposed into an orthogonal sum of the spaces 
Hq. Hence, any F G L‘^{Q) has a IFzener chaos expansion 

00 

f = YW,). 

5=0 

(11) 

where /o = E[F] and the fq G 77®“^, g > 1 are uniquely determined by F ( 
1 .1.2 of [H]). 

'see Theorem 
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The purpose of the above discussion is to provide sufficient background to use the 
Fourth Moment Theorem. This theorem, hrst published in 2005, has inspired an exten¬ 
sive body of literature, and provided solution techniques to a new class of problems. This 
hrst version of the theorem was proved in [16]. Since then, other equivalent conditions 
have been added [lailS]. A key advantage of this theorem is that, unlike the method 
of moments, it is sufficient to check the convergence of the moments of up to order four. 

Theorem 2.1 (Fourth Moment Theorem). Fix an integer q > 2. For integers n > 1, 
let Fn = Iqifn) 0 Sequence of random variables belonging to the qth Wiener chaos of 
X, so that fn G Assume that E [Fff\ —)■ cr^ > 0 as n —)■ oo. Then the following 

are equivalent: 

(a) As n ^ oo, Fn converges in distribution to N A/'(0, cr^). 

(b) lim„^oo E [F^] = 3a" = E [iV"]. 

(c) For each integer I < r < q - 1, hm„_^oo ||/n <8)r /n||^® 2 ( 9 -r) = 0. 

We have this multidimensional extension due to Peccati and Tudor: 

Theorem 2.2 ([H]). For an integer k > 1, let Fn = {Ff,...,F^) be a sequence of 
random vectors, and let 1 < qi <■■■< qk < oo be integers such that for each 
j = Fi = Iqj(fi) for some kernel ff G "H®®. Moreover, assume that each 

E [{Fl)^] —y aj >0 as n tends to infinity, and that lim^^ooE = 0 for all j ^ 1. 

Then the following are equivalent: 

(a) Fn converges in distribution as n ^ oo to S), where S is the k x k diagonal 
matrix with entries 

(b) For each j = 1,..., k, Ff, converges in distribution as n ^ oo to Af{0, a|). 

For a general sequence of square integrable random variables, the convergence to the 
normal distribution can be deduced from the Wiener chaos expansion. Using Theorem 
12.21 one can show that if every projection on the Wiener chaos satisfies the hypotheses 
of the Fourth Moment Theorem, then their limits are independent and a central limit 
theorem holds for the global sequence. This phenomenon can be described as a chaotic 
central limit theorem. 

Theorem 2.3 ([H])- Let {Fn} be a sequence in L^(r2) such that E [Fn] = 0 for all n. 
Write each Fn in the form 

OO 

Fn = Iq{fn,q)i 

q=l 

and suppose that the following conditions hold: 

(a) For each q>l, lim^^oog! ||/n,q ||^®9 = for some a^ > 0. 

(b) a^ = i:7=i<<^- 

(c) For each q>2 and r = l,...,q-l, lim„^oo \\fn,q /n,g||^® 2 ( 9 -r) = 0. 
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(d) lim^^ooSup„>i \\fn,q\\lim) = 0- 

Then as n tends to infinity, —> ■A/'(0, a^). 

In the following sections, the symbol C denotes a generic positive constant, which 
may change from line to line. The value of C might depend on T and the properties of 
the process X. 

3 Central limit theorem for variations of a self-similar Gaussian 
process 

3.1 Defining characteristics of the process 

Let X = {Xt,t > 0} denote a centered self-similar Gaussian process with self-similarity 
parameter fi G (0,1). We introduce the following conditions on the function dehned 
in ([T]), where a G (0,2/3]: 

(H.l) (j) has the form (j){x) = — X{x — 1)'^ + 'ijj{x), where A > 0, 'ijj{x) is twice-differentiable 
on an open set containing [1, cx)), and there is a constant C > 0 such that for any 
X G (1, cxo) 

(i) \'fi'{x)\ < Cx°'~^; 

(ii) |'^"(x)| < Cx~^{x — 1)““^; and 

(iii) = /3'0(1)) when a > 1. 

The lemma below shows that a satishes E [(W+g — Xtfi] ~ s“, and we call this value 
the increment exponent. 

Lemma 3.1. Under (H.l) for 0 < s < t we have 

E [(W+. - Xtf] = + gfit, s), 

where | 5 fi(f, s)| < CsU^~^ if a <1 and | 5 fi(f, s)| < Cs‘^U^~‘^ if a>l, for some constant 
C. 

Proof. We can write 

E [(X,+, - A',)2] = (( + s)'».^(l) + ("-’.^(1) - (l + 

=^( 1 ) ((«+sr - n + 2^= (^( 1 ) - ^ (i +/)) 

= .^(1) (((+ sff - + 21^0 (v.(i) + ^ (j)“ - v* (i + /)) 

= 2Af2^-“s“ -f V'(l) {{t + - t^^) - 2^^ / fjfiy) dy 

= 2At2/5-"s" + ^i(f,s). 
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If a < 1, it follows from the Mean Value Theorem and the estimate on |'0^(2^)| that 
s)| < . If a > 1 we use property (hi) to get 

s) = - ^|J'iy)]dy, 

which implies that \gi{t, s)| < due to the estimate on |' 0 "(a;)|. □ 

Notice that Lemma [3.11 implies that for 0 < s < t, 

( 12 ) 

for any e > 0 such that 1 — a — e> 0 ifQ!<l and 2 — a — £> 0 ifl<a< 2 . 

Lemma 3.2. Assume condition (H.l). Then, 

(a) For 0 < 2s < t, we have 

E [(Xt+, - X,){X, - X,_,)] = (2“ - 2)At2^-“s“ + g 2 {t, s), 
where \g 2 it, s)| < Cs‘^(t — + Cs“+^(t — 

(b) For 0 < 2s < I < r < t — 2s, 

E[(Xi-Xi_,)(X,-X,_,)] 

= A(r - s)^^“" [(t - r - s)“ + (t - r + s)" - 2{t - r)"] + g 3 {r, t, s), 
where \g‘i{r, t, s)| < Cs^(r — s)^^“““^(t — r — s)"“^ + Cs^(r — s)^^“^. 

Proof. For (a), we can write E [(V^+s — Xt){Xt — V'i-s)] as 
E [Xt+sXt - X^ - Xt+sXt-s + XtXt-s] 

= i^'>(0(i + £)-^(i))-(i-.r(^(i + ^)-^(i + ^)) 

= _ s) 2 / 3 -a (A(2s)“ - As“) 

+ ^''(l + y) dy-{t- sf^ V'' (^1 + dy 

= A(2“ - 2)t^X»s^ + (1 - 2")As" - (t - sfX») 

+ -{t- sff^) + y) dy 

+ (t-s)^^^ ^|J'{l + y)-Ip'(^1 +j^+y^ dy 

= (2“-2)As“t2/3-« + ^ 2 (t,s), 


+ y]dy 








where, given conditions on ^|J and its derivatives, we have that 
\ 92 {t,s)\ < Cs\t - + Cs^+\t - 

Next, for (b). 


= 




= Ar2^-“ ((t-r- s)“ -(t- r)“) - A(r - ((t - r)“ - (t - r + s)“) 


+ r 


2/3 




t — s 


+ y] dy-{r- s) 


2/3 




= A(r — sf^ “ ((t — r — s)“ + (t — r + s)" — 2{t — r)") 

- A - (r - ((t - r)“ - (t - r - s)") 


+ r 


2/3 


V'' 


t — s 


+ 3/ dy-{r- s) 


2/3 




t — S 


r — s 


t — s 


r — s 


+ y] dy 


+ y] dy 


where 


= A(r - sf^ “ ((t - r - s)“ + (t - r + s)" - 2(t - r)") + g-i{r, t, s), 

gsir, t,s) = X ((r - {{t - r)“ - (t - r - s)") 

t — s 


+ - (r - 

t — s 




r — s 


+ y]dy 


,2/3 


_ r-2/3 




r — s 


+ y] 


, f t — s 


+ y 


dy 


f ( -—- + 3/1 dy. 


Using the Mean Value Theorem, the fact that | < r—s < r and bounds on the derivatives 
of I//)I, we have that, for a constant C, 


|/? 3 (U t, s)| < Cs^r - - r - s)""' + Cs\r - 


t — s 


r — s 


0—1 


r — s J \r[r — s) J 



+ CsHr-s) 


If a < 1, then {t — s)“ ^ < (t — r — s)“ ^ and {t — r)“ ^ < {t — r — s)" Therefore, we 
have the bound 


|// 3 (r, t, 5) I < Cs\r - - r - 
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In the case a > 1, then r > t/3 implies {t — s)/{r — s) < 5. Using this inequality and 
t — r < 2{t — r — s), we obtain 

Mr, t,s)\< Cs\r - - r - + Cs\r - 

Then, the proof of part (b) is complete. 

□ 

Example 3.3. Let > 0} denote a fractional Brownian motion with Hurst 

parameter H. Then condition (H.l) is satished with a = 2(5 = 2H, A = |. In this case, 
we obtain 

E [(B«, - B«r] = 

E [(B,+. - - B»,)] = i( 2 "« - 2)5^" 

E [(B,» - - B",)] = i ((( - r + sf" + (( - r - - 2(t - rf «), 

which means that = 5^2 = 5'3 = 0 in Lemmas 3.1 and 3.2. This means that, in the 
general case, we can think of X as a process that is similar to the fBm, but with an 
additional, lower-order correction term on the covariance. In Section 4 we give examples 
where the terms pi, i = 1,2, 3, are nonzero. 


We will make use of the following additional condition on the behavior of the hrst to 
derivatives of 0 at inhnity, which cannot be deduced from condition (H.l). 

(H.2) There are constants C > D and 1 < v <2 such that for all x>2, 


(i) 
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| 0 '(a;)| < 
|0"(x)| < 


C(x- 1)-^ 
C{x - l)“-2 

’u(x-l)-^-i 

C{x - l)“-3 


if a < 1 
if a > 1, 

if a < 1 
if a > 1 


3.2 Central limit theorem 

We are now ready to state the main result. 

Theorem 3.4. Suppose a self-similar Gaussian process (X, 0) satisfies (H.l) and (H.2) 
above. For n >2, consider the stochastic process defined in that is, 

\nt\-l 

Fn(t) = ^ /(^j,n), \nt\ > 1, 

i=o 

and Fn{t) = 0 if \nt\ < 1 . ITe assume that f G L^(M, 7 ) has the expansion ([^ with 
Hermite rank d > 2. Then, if a < 2 — A, the finite dimensional distributions of the 
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processes {Fn,n > 2} converge in law to those of a Brownian motion with scaling given 
by where 

ctJ = 2“V! + ^1“ + 1"^ - ir - 2|m|“)''. (13) 


3.3 Proof of Theorem 13.41 

We show the convergence of the hnite dimensional distributions using Theorem 12.21 and 
Theorem 12.31 For an integer p > 2, choose times 0 < ti < • • • < tp < oo, and dehne 
to = 0. For i = 1,... ,p, dehne Gn{ti) = Fn{ti) — We want to show that each 

Gnifi) converges in law to A/'(0,(j^(ti — ti-i)) applying Theorem 12.31 Without loss of 
generality, it is enough to prove this for Gn{ti) = Fn{ti). 

The projection of Fn{ti) on the Wiener chaos of order q is 


[ntij-l 

i=o 

By (HQD, we can write Fn^q in terms of the stochastic integral operator Iq 


IntiJ-l 



where we use the notation 

ij,n = ||^^i||L2(0)- 

71 

The symbol dj/n denotes the indicator function of the interval that is di = 

l[i. i±li- To verify the conditions of Theorem 12.31 we adopt the following Hilbert space 
notation. The indicator function l[o^t] is an element of the Hilbert space dehned as 
the closure of the set of step functions with respect to the inner product 

(l[o,s], l[o,t])y^ = E , s, f > 0. 

With this representation, we can write Xn,(?(H) = Iq{fn,q), where 


[ntij-l 




i,q — fn,q{tl) — ^ ^ 

^ j=0 




(14) 


where hj^n = B is apparent that /„ G 

Now we proceed to verify the conditions of Theorem 12.31 

Proof of condition (a): We want to show that for any g > d, E [Fn^q{tiY] converges to 
(Tqti, where is given by (IT^ . This is the contents of Lemma [5.21 whose proof is given 
in Section 5. 
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Proof of condition (b): This is obvious by definition. 

Proof of condition (c): We wish to show that for each r = 1, 2,..., g — 1, 

lini \\fn,q^rfn,q\\f.^2(q-r) 0. (fb) 

n—j-oo ^ 

By ([8]) we have for each r = 1,..., g — 1 


^2 KlJ -1 

fn,q fn,q ^ ^ ^fc,n) 

j,k=0 


I {hfr 


k.n 


0 . 


which is an element of It follows that 


^4 KlJ -1 

\\fn,q^r fn,q\\f^(Si2(q-r) — ^ ^ {hl^ni hm^n)^ {hk^W) hm,n)^ 

j,k,£,m=0 

( 16 ) 

We proceed in a manner similar to the proof of the convergence (1.3) in [S]. Fix an 
integer M > 1. First we decompose the set of multi-indexes D = [0, [ntij — 1]'^ into 
D = Di^m U where 


fc, £, m 

)eD-.\3-k\<M,\l- 

m\ < 

M, \j-e\<M}, 


D2,m — D^,m U Di,M U Dq^m, 


Dz,m 

= {{j,kj,m) e D :\j - 

■k\ > 

M}, 

Da^m 

= {U,k,i,m) e D :\e- 

m\ > 

M}, 


= {{j,k,i,m) e D :\j - 

-i\> 

M}. 


and 


Taking into account that = 1 and using Cauchy-Schwarz inequality, it follows 

that 






q—r 

Sj I - 


< 


CM^ 


{j,k,e,m)GDi^M 


n 


which converges to zero as n tends to infinity, for any fixed M. It suffices to handle the 
sum over one of the sets Di^M for * = 3,4, 5. The analysis is the same for each of them 
and we consider only the case i = 3. Set 


^n,M — 


4 E 


^ 


'm^n/ 


\q-r\ 
/i 
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By Holder’s inequality, we can write 


^ 2 1 ^ ^ I I | (^£,n) 

{j,k,i,m}eD3 

,M 


1-2 


X 


< 


^ ^ I (^_7,n5 I hm,n) ^ 

0 >,£,m)GD 3 

,M 


^ ^ NiJ-i 

^ ^ I (,hj,ny hk,n} 


n 


\ 


j,k = 0 


[nliJ -1 


2-22 

<J 




n 


E K'* £,n) 


£,m=0 


\j-k\>M 

Using the same arguments as in the proof of Lemma [5.21 one can show that the lini sup 
as n tends to inhnity of the above expression is bounded by 


2 -, 

mgZ,|m|>M 


|m + 1|" + |m — 1| 


X 


2 ||m + 1|“ + |m — 1| 


2 |m 


2 m 


mGZ 



which converges to zero as M tends to inhnity. This completes the proof of property 
(c). 

Proof of condition (d): Condition (d) follows from the proof of Lemma 15.21 In fact, 
we know that the series c^q\ is convergent, and it suffices to take into account the 
estimates (|2T|) . (i26|) and Remark 5.3 after the proof of Lemma [5.21 

With conditions (a) - (d) of Theorem 12.31 satished. it follows that for i = 

Gn{ti) converges in law to a Gaussian random variable with mean zero and variance 
given by — U_i). We next want to show that lim„_,.oo E = 0 when 

ti 7 ^ tj. Without loss of generality, it is enough to show that 


sup 

q>d 


[ntij-l [nt 2 j-l 

-EE & 


-q^-q 
j,n 


k=0 j = [ntiJ 


/ 


\ 9 

IE 

AX^AXk 

) 

V 

n n . 

/ 


tends to zero as n tends to inhnity. Let = {nt 2 ) s, we can decompose the above sum 
into 


^ E & 


-q ^-q 

'j,n 


{j,k)eDi^, 


E 


AX, AXi 


n n 


(j,A:)eD 2 ,n 


E 


AX, AXi 


n n 


where 


Di,n = {(j, k) e Dn: \j - k\ < M„}, 


D2,n = {(j, k) e Dn: \j - k\> Mn}, 
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and 

Dn = {{j, k) :0<k < lnti\ - 1, [ntil < J < ^^ 2 ] - !}• 

By Lemma [3.II and Cauchy-Schwarz, the first sum is bounded by Cn~^. For the second 
sum, using arguments from Step 4 and Step 5 in the proof of Lemma [5.21 (see fl25]) ). the 
lim sup as n tends to infinity is bounded by 

00 

C lim V ((m + 1)“ + (m - 1)“ - 2m^y < C lim = 0, 

n—>-oo ^ ^ n—)-oo 

m=Mn 

where the limit follows since = Cn^ and q(a — 2) -|- 1 < 0. 

Based on the above. Theorem [2]2] implies that the vector sequence (G„(fi),..., Gn(tp)) 
converges in distribution to A/'(0, S) as n tends to infinity, where S is the diagonal matrix 
with entries — tj-i), i = 1,... ,p. Taking Fniti) = X]L=i Gnitk), we have that 

(F„(ti),..., F^itp)) A (F(F), • • •, F{tp)) (17) 

as n ^ 00 , where each F{ti) ~ A/'(0, This completes the proof of Theorem 13.41 

3.4 Theorem 3.4 examples 

In this section we consider some particular cases where the function / can be repre¬ 
sented as a single Hermite polynomial, or a linear combination of finitely many Hermite 
polynomials. 

Example 3.5. For integer p > 1, let f[x) = — E where Z has a A/'(0,1) distri¬ 

bution. The Stroock formula [12] for a Hermite expansion gives: 

2P 

/(Z) = E [/(Z)] + -E [/W(Z)] Hyz) 

q=l 

= -^E \Z^P-y HJZ). 

“ q\{2p — q — 1)\ ^ 

Since 2p — 1 is odd, / has Hermite rank 2, and Theorem 13.41 can be applied. In this 

case, Fn converges in law to a Gaussian random variable with variance 

where 

^ = (2p)K2p - <1 - i)n 

^ q\{2p — q — l)\ g!(2p —g —1)! 

for even q = 2,, 2p, and Cq = 0 for all odd integers and q greater than 2p. 

For the odd integer case, let f{x) = — E [iZp^+^j. Again, since E [f'{Z)] = 0, 

the Hermite rank is 2 and we can apply Theorem 13.41 In this case, 
where 

^ (2p+ 1)! n |2p_g+n ^ (2p + l)!(2 p-g)!! 

^ q\{2p — q)\ ^ ^ q\{2p — q)\ 

for even g = 2 ,... 2p, and Cg = 0 for all odd integers g and all g greater than 2p. 
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Example 3.6. Suppose f = Hg for a single Hermite polynomial of order q > 2. In this 
case, one can show that for fixed t, Fn = Fn{t) converges in total variation. Let N 
denote a random variable with the J\f{0,tag) distribution. In [12] it is proved that 


d-Tvi^n, N) < 






where cItv denotes total variation distance and DFn is the Malliavin derivative of Fn 
[I31[I1]. We can write Fn = Iq{fn,q), where fn,q is given by fl20|) with ti = t. From [T^ 
Lemma 5.2.4], we have 


Var 


DFn 


q-l 

r=l 


t^t\ 


{2q-2r)\ 





2 


where fn,q fn,q denotes the symmetrization of fn,q 0r fn,q- Therefore, by condition 
(c) in the proof of Theorem 13.41 and using the identity 

2 2 

^ — \\fn,q’^rfn,q\\sj2(.q-r) 

for each r = 1 ,..., g — 1, we obtain 




n^q 




n,q 


lim dTv{Dn,N) = 0. 

n—)-oo 

On the other hand, in this example we can show a functional central limit theorem. 
Indeed, by (|27|) in the proof of Lemma [5.21 and Remark 5.3 following, we can show that 
for sufficiently large n, 

E [{Fniu) - Fn{U.i)f] < C . 

Moreover, using the fact that all the p-norms are equivalent on a fixed Wiener chaos, 
for arbitrary 0 < Si < s < S2 < T we deduce that for n large enough 

E [{Fn{s) - Fn{si)f {Fn{s 2 ) - 

< (E [{Fn{s) - Fn{s,))^]y (E [{Fn{s2) - 

This implies that the laws of the processes {Fn,n > 2} are tight in the Skorohod space 
D[0, 00) (see Billingsley [H Theorem 13.5]). As a consequence, from Theorem 13.41 we 
deduce that the laws of Fn converge in the topology of D[0, 00) to a Brownian motion 
with scaling a^. 
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4 Examples of suitable processes 

4.1 Bifractional Brownian motion 

Bifractional Brownian motion is a generalization of fBm, first introdnced by Hondre 
and Villa in [H]. It is defined as a centered, Ganssian process B = ,t > 0} with 

covariance 

E [BA] = i [(s2» + e")‘‘ -\t- . 

where H G (0,1) and K G (0,1]. Note that if it' = 1, then B is an ordinary fBm. The 
reader may refer to [ini[IB] for a discussion of properties. 

The covariance can be expressed in terms of 0 with {3 = HK and 

^ [(1 + • 

As stated in Section 1, it is well known that for ordinary fBm, Fn{t) converges in 
distribution for all iJ < 1 — Since the case it' = 1 is well known, we will assume 
below that i^ < 1. We now verify the properties (H.l) - (H.2). 

For (H.l), we write 

- ir'"+^(1+ 

which means we have A = 2~^ and a = 2{3 = 2HK. Then = 2“^(1 + , and 

the bounds on l"^'! and l"^"! follow. Moreover, = HK = /3'0(1). 

For (H.2), 

(t)'{x) = 2^-^HK {x^^-\l + -(x- i)2-f^^-i) . 

We can write 

(1 + + (A' - 1) f\x^^ + yf-^dy, 

Jo 

so that 

0'(x) = 2^-^HK - (x - + (iF - + y)^-'^dy^ . 

Hence, we can write |0'(x)| < C{x — lY^^~‘^ + C{x — lY^^~‘^^~^. If a < 1, then we take 
V = min {1 + 2H — 2HK, 2 — 2HK} > 1, and if a > 1 then 2 — 2HK <2 + 2H — 2HK 
implies |(/>Xx)| < C{x — i)2iii^-2_ Continuing, we have 

4)'\x) = 2^-^HK{2HK - 1) - (x - l)2-f^^-2) 

+ 2‘^-^H‘^K{K - l)x^^-‘^ [ (x^^ + y)^-^dy 

Jo 

+ 2‘^-^H^K{K-1){K-2)x^^-^ f {x'^^ + yf-^dy, 

Jo 
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so it follows that 

W’{X)\ < C{X - + C{X - lfHK-2H-2 < ^ 

where 

min{2 + 2H — 2HK, 3 — 2HK} = z/ + l ifQ!<l 
3 - 2HK if a > 1 ■ 

4.2 Subfractional Brownian motion 

Another variant on the fBm is the process known as sub-fractional Brownian motion 
(sfBm). This is a centered Gaussian process {St,t > 0}, with covariance defined by: 

fl(s, t) = - i [(s + f)“ + Is - ir"] , (18) 

with real parameter H G (0,1). Some properties of the sfBm are discussed in [2113]. Note 
that TT = 1/2 is a standard Brownian motion, and also note the similarity of R{s,t) to 
the covariance of fBm. Indeed, in |1] it is shown that sfBm may be decomposed into an 
fBm with Hurst parameter H and another centered Gaussian process. 

Let S = {Stit > 0} denote a sub-fractional Brownian motion with 0 < iL < 1. In 
this case we have A = i, a = 2/3 = 2H and 

(/)(x) = 1 -t- ^ ((x (x - 1)^^) = -A(x - 1)^^ -F 'i/(x), 

where 'ip{x) = 1 -f x^-^ — |(x -|- l)^^. Glearly, \'iIj'{x)\ < Cx‘^^~^ and \'il)"{x)\ < Cx‘^^~^. 
Moreover, = H{2 — 2‘^^~^) = (3'tp{l). 

To check (H.2), we have 

|(/)'(x)| = H\{x + - 2x‘^^-^ + (x - < C{x - 

so (H.2)(i) is satisfied for x > 2, and (H.2)(ii) can be shown by a second derivative. 

4.3 A Gaussian process introduced by J. Swanson 

We consider the centered Gaussian process Y = {Yt,t > 0} with covariance given by 

E [YsYt] = Vsf sin“^ ’ 

Then Y can be characterized as a self-similar Gaussian process with (3 = 1/2 and 

(/)(x) = v^sin 

This process was studied by Jason Swanson in a 2007 paper [20], and it arises as the limit 
of normalized empirical quantiles of a system of independent Brownian motions. The 
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properties of this process were also considered in Unlike the fBni family processes 
in Sections 4.1 and 4.2, this process is an example of the case a < 2/3. 

This process satisfies (H.l), with /3 = 1/2 and a = 1/2. We can write this as 

(j){x) = —\/x — 1 + ^-y/xsin”^ ~ ^ ’ 


which gives A = 1 and 


ipix) = \/xsin 


-1 


+ \/x - 1. 


Then 


W{x)\ = 


X 2 


sm 


T — 1 


Vx -1 


< Cx 2 . 


And for the second derivative, 


r{x) = 


_ 3 
X 2 


sm 


-1 


v^- 1 


_1 
X 2 


\/x-I J ^s/x - 1 V\/^(\/^ + 1) / ’ 
so |'?/"(x)| < Cx~^ ^1 + (x — l)“^j < Cx~^[x — 1)“^. 

To check condition (H.2), we return to the original expression (j){x) = y^sin“^ (^x~^ 
We have 


[X = 


2VT 

1 


sm 


^/x — 1 
1 


+ 


sm 


-1 


2^/x \y/x y/x -IJ 2yfx 
From a Taylor expansion of sinf, we have for 0 < f < 1 

for some 0 < h < f, and it follows that 

■if 

We then set t = x~^ , and use a second Taylor expansion on sin”^, to conclude that for 

X > 2, 


sm 


-1 


< 


C 


Hence, we have 


|0'(x)| < 


2^ 


X y/x-1 


+ 


Xy/x — 1 

c 


X2 y/x — 1 


^^(x-l)-^ 
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so we have v = 2 since a = | < 1. 
Similarly, 



Hence, |0"(x)| < C{x — 1) ^ = C{x 


\/x — 1 
1 


—7 a; 2 

4 


x\/x — 1 
1 


y/x — ij 

X > 2. 


X 


(x — 1) 2^ 


4.4 Two smooth processes with a < 1 


For 0 < a < 1, we consider the centered Gaussian processes Zi{t), Z 2 {t), with covari¬ 
ances given by: 

E [Zi{s)Zi{t)] = r(l - a) ((s + - max(s, t)“) 

E [Z 2 is)Z 2 it)] = r(l - a) (s“ + - (s + t)“). 

These processes are discussed in a recent paper by Durieu and Wang j6], where it is 
shown that the process Z = Zi + Z2 (where Zi, Z2 are independent) is the limit in law 
of a discrete process studied by Karlin. The process Zi is new, but the process Z2, with 
a different scaling constant, was first described in Lei and Nualart nni. 

For Zi, we can write E [Zi{s)Zi(t)] = s°‘(j)i{t/s), where 01 (x) = F(l—a) ((x -f 1)“ — x“). 
We can also write 0i in the form 


0i(x) = -F(l - a)(x - 1)“ F(1 - a) ((x 1)“ (x - 1)“ - x“), 

so that A = F(1 — a) and 01 (x) = F(1 — a) ((x -1- 1)" -|- (x — 1)“ — x“). Note that by a 
Taylor expansion for x > 1, 


(x — 1)“ = x“ — ox" 


0-1 I 1 ) «-2 


-X 


+ O(x“-0, 


(19) 


so that we have |0((x)| < Gx“ ^ and |02Xx)| < Gx“ and Zi satishes (H.l). For (H.2), 
note that we can write 

01 (x) = aF(l — a) f {x + u)°‘~^du, 


so for X > 2 

|0'^(x)| = a|a — 1|F(1 — a) j {x + u)°‘~'^du < 2x°‘~‘^, 

Jo 

which satisfies (H.2)(i) with u = 2 — a > 1; and it can also be seen that |0"(x)| < 
4x«-3 _ 2x“^“^, so that (H.2)(ii) is satisfied. 


For Z 2 we have 


02(x) = F(1 - a)(l -f- x“ - (x 1)“) 

= -F(l -a){x- 1)“ + F(1 -a){l + x^ + {x- 1)" - (x + 1)") 
= -F(l - a){x - 1)" 02(a;), 
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where again we take A = r(l — a). By a computation similar to that for Zi above, it 
can be seen that "02 satisfies (i) and (ii) of (H.l). The computations for (H.2) are also 
similar to those for Zi above. We write 

= T{1 — a) — aV{1 — a) f {x + u)^~^du, 

Jo 

so that (H.2) conditions (i) and (ii) are satished with u = 2 — a > 1. 

Remark 4.1. Both processes satisfy E[Zj(f)^] = Cf", so the increment exponent at 0 is 
a. On the other hand, for t > 0, E [{Zi{t + s) — Zj(f))^] qs for a constant q, and the 
increment exponent for f > 0 is 1. The renormalization in Theorem 13.41 is given by a, 
which is also the self-similarity parameter. 


5 Some technical lemmas 


We begin with a technical lemma that gives upper bounds on certain covariance terms. 

Lemma 5.1. Let n > 6 be an integer, and let j,k > 1 be integers satisfying 3k < j. 
Then under (H.2), there is a constant C > 0 such that 

— k)~'' if a < 1 

Cn~‘^^k'^^~°‘(j — if ct > 1 ' 

where the exponent 1 < u < 2 is defined in (H.2). 

Proof. We have 

'j + 1' 


E 

AX^AXk 

<i 


n n _ 

^ 1 


E 


AX,AXk 


= n-^^ 


{k + lf^ ( 0 


k + 1 


' 3 + 1 


3 

k+1 

I 
k 


3 + 1 
k + 1 


k + 1 




3 

k + 1 


J + 1 
k 


+ ' 


k + 1^ 

Condition j > 3k and (H.2) imply that there exists a constant C > 0 such that for each 
i±l-\ 

^ Lfc+i’ fc+iJ 

'Ck^{j-k)- 

ce-^{j - ky-^ 


\ct>'{x)\< 


if a < 1 
if a > 1 


and for each x e , 




-u-l 

0—3 


if a < 1 
if a > 1 


Then the desired estimate follows easily from the Mean Value Theorem. 


□ 
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Lemma 5.2. Let t > 0. Under above definitions with a < 2 — ^ 

lim E [Fl it)] = aj, 

n—^OO L J >1 

where is given by oa. 

Proof. The proof will be done in several steps. 

Step 1. It follows from (|9]) that 


2 \nt\-l 

E [f;h(<)] = 1 E E Iff, (r,,n) ff, (n,„)i 

j,k=0 


ff E ii^Eiiz’(«) 

j,k=0 




AX, AX* 


1 \ <1 


n n 


Fix 7 G (0,1/2) and decompose the above double sum into two terms, that is, E [F^g] = 
Ai,n + ^2,n, where 


^ II II 1,2(0) 

j,k&Di 


||AXi| 


L2(o) 


E 


AX, AX* 


1 \ <1 


n n 


and Di = {j : 0 < j < n"^ A i[nt\ — 1)} and D 2 = {j : n"^ < j < [nt\ — 1}. 

The term Ai^n can be bounded, using Cauchy-Schwarz inequality, by and 

it converges to zero as n tends to infinity. So, it suffices to consider the term A2,n- We 
recall the notation = ll^Wj ||i2(Q). 

n 

Step 2. From we can write for j G D 2 



[l + V,,n], 


where 


\Vj,n 


2X-’ 


9i 



Consider the following decomposition; 




( 20 ) 

( 21 ) 


_ 


'“‘Z IE 


AX, AX* 


= 


j,kGD2 

[nt\ - 1 - q\Cg 


+ 


n 


n 


E ZAZ (e 


j,kGD2,\k-j\=l 


AX, AX* 


n 


E zii 


-QC-1 

k.n 


j ,k^D2,\k—j\>2 


E 


AX, AX* 


The first term clearly converges to q\c^gt. We denote the second and third term by i?2,n 
and i?3,n, respectively. 
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Step 3. Let us consider the term -B 2 ,n- Using Lemma [3.2f ai we can write 

1 1 


E 


AX,AXj-i 


= A(2" - + g2 


n n 


Multiplying this expression by using fl20|) yields 




AXj AXj-i 


— [1 + rjj^n] 2 [1 + Vj-l,n] ^ 


X (2“-^ - 1) 


i -1 


i(2/3-a) 


+ R 


'J,n I 5 


( 22 ) 


where 


n n 


= {2X)-V^{j{j - l))-5l- 

Applying Lemma [3.2r ai and assuming > 2, this term can be bounded as follows. 

<^-(2/3-“+^), (23) 

where 5 = min(2 — 2/3 ,1 + a — 2/3). We claim that there exist p G (0,1) such that for n 
large enough and for all j E D 2 , 


C',nC-\,nE 


AX,' AX,-i 


< p < 1. 


(24) 


This follows from the estimates fl2T]) . fl2^ and the fact that |2" ^ —1| < 1 and 2/3—a+5 > 
0. Finally, from the expression fl22l) and the estimates fl?Il) and fl23l) it follows that 

lim B 2 ,n = q\cl2 lim - [1 + Pyn ]”2 [1 + 

^ fi^oo Tl ' ^ 


ij-ier >2 


X 1^(2"-^ - 1) 
= g!cj 2 ( 2 “-^ - l)''f. 


i - 1 


i(2/3-«) 


Rj,n 


(25) 


Step 4- Let us consider the term B^ n- First we show that the terms with k < [j/3j or 
j < L^/3j do not contribute to the limit. That is, we claim that the following expression 
converges to zero as n tends to inhnity 


[ntj-l [3J 

= E T.&. 


-q^-q 
k,n 


j =3 k=l 


E 


AX,' AX* 


n n 


To do this, we consider two cases. 

Case 1. When a < 1, Lemma [5.11 gives 


\nt\-l [§J 

Cn<Cn-^ 

i=3 k=l 




















22 


Note that 1 < /c < j/S implies (j — k) '^ < Cj so it follows that 

„ \nt\-l [Ij \nt\-l 

^ i=3 fc=l ^ j=3 

which converges to zero as n tends to infinity because q{a — 2) + 1 < 0 since a < 2 — ^ 
and q > d. 

Case 2: Assume 1 < a <2 — For this case, by the other part of Lemma [5.11 


C 


\nt\-l [§J 


Y1 


q{a-2) 


j=3 k=l 


Again, since 1 < A; < j/3, we can say that {j — k)^ ^ < Cj°' where C* is a constant 
that does not depend on j or k. Thus, we have 


C 


[ntj — 1 


q(a-2)+l 


Cn<-J2l 

j=3 

which converges to zero as n tends to infinity because q{a — 2) + 1 < 0 since a < 2 — ^ 
and q > d. 

Step 5. Finally, it remains to study the following term 

2qK 


Dn = 


n 


E 


(E 


AXjAX, 


n n 


j,keD2,[j/3\<k<j-2 


We have by Lemma [3.2l bL 


E 


AXiAX, 


n n 


= Xn '^Aj^k + 93 


k+l j+l 1 


n n n 


where 


Aj,k = U - + 1|“ + |j - fc - 1|“ - 2\j - 

Multiplying this expression by using fl2U]) yields 

'k + l j + l 1 


- 1^-1 




•j,n^k,n 


k + 


93 


n n n 


— [1 + 2 [1 + rjk,ri\ ^ 

X ( 2 -^ {k/3Y^^^-^X,^k + Rpk,n) , 


Rj,k,n = ^n‘^^{kj) 2 ^^^ “^ 5(3 


k+l j+l 


n 


n 



where 
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We can write 


r=l ^ ^ 


From the estimate for l^fal from Lemma [3.2r bh we have 

\Rj,kA < C'j-^(2/3-a)^i(2/3-a)-l^^. _ ^ ^ ^^.-l(2/3-a)^i(2/3-a)+(a-2) _ 


It follows that we can write 


r=l 

Using that k/j < 1, \AjA E (FIj —/c —1|““^ for |j —/c| > 2, /c“'’ < C^j~'^ < C"’(j —fc —1)“'’, 
and similarly that /c““^ < C{j — k — 1)"“^, we obtain 

for some constant C > 1. This implies 


2(7'c^ ^ ^ 2rj'r^ 

—^ ^ [^ + Vj,n]~H^ + Vk,n]~^^2,j,k,n\ < —^C'^ Y1 


2 L’T’d-i 


n 


n 


j=[n'Tj fc=[j/3j+l i=F'^J 

which converges to zero as n tends to inhnity. Moreover, there exists p G (0,1) snch that 

\nt\ —1 


1 


n 


c- E i 

j = [nT'J 


q{a-2)+l 


< P 


(26) 


for n large enongh. 
Therefore, 


lim Dn = q\cl lim - V 2-'?(fc/j) 

j,keD2,lj/3i<k<j-2 

Since we know from Step 1 and Step 4 that the terms with j, k E Di and k < [j/3j do 
not contribnte to the limit, we can add these terms and write 

\nt\-l j-2 

lim Dn = q\A lim - V V 

n—>-oo ^ n^oo ft ' ^ ^ ^ 

i=3 k=l 

We nse the change of index m = j — k to write 


[nij—1 [nt\—k—l 


lim Dn = q\cl lim — 2 


k=l m=2 


k + m 


f(2/3-a) 


{{m + 1)“ + (m — 1)“ — 2m")^ , 

( 27 ) 
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which leads to 


lim Dn = 2^-^q\clt V ((m + 1)“ + (m - 1)“ - 2m")''. 

1.^00 ^ 


(28) 


Thus, the limit (fT3|) follows from results (|2^ and (|28|) . and Lemma [5.21 is proved. 

□ 


Remark 5.3. Notice that the right-hand side of (127p is bounded by 


2q\c]t Y, 2-“ ((m + 1)“ + (m - 1)" - 2™“)" 

m=2 


where sup^>2Pm < L 
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